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Abstract: We determine the equilibria of a rigid loop in the plane, subject to the constraints 
of fixed length and fixed enclosed area. Rigidity is characterized by an energy functional quadratic in 
the curvature of the loop. We find that the area constraint gives rise to equilibria with remarkable 
geometrical properties: not only can the Euler-Lagrange equation be integrated to provide a quadrature 
for the curvature but, in addition, the embedding itself can be expressed as a local function of the 
curvature. The configuration space is shown to be essentially one-dimensional, with surprisingly rich 
structure. Distinct branches of integer-indexed equilibria exhibit self-intersections and bifurcations — a 
gallery of plots is provided to highlight these findings. Perturbations connecting equilibria are shown to 
satisfy a first order ODE which is readily solved. We also obtain analytical expressions for the energy as 
a function of the area in some limiting regimes. 
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1 Introduction 

Our object of study in this paper is a deflated rigid loop {elastica hypoarealis) . To explain this, consider 
a closed loop in a plane. The loop is made of some elastic material and, if taken thin enough, the 
contribution of its longitudinal deformations to its elastic energy can be ignored. This gives rise to a 
hamiltonian density proportional to the square of its curvature — hence elastica. It is clear that fixing 
only the length L of the loop will not lead to any surprises, the only possible equilibrium being a circle. 
What will happen though if we, in addition, fix its area A to be less than that of the circle? (hence 
hypoarealis) . 

This particular hamiltonian (without the constraints) has a distinguished history, dating back at least 
to the BernouUis. It makes a reappearance in Euler's inspired analysis of a mundane but subtle problem 
in mechanical engineering: the buckling of a loaded beam |2^, ^. More recently, variants of this problem 
have attracted the attention of mathematicians. In particular, the problem of determining the curves of 
constant length which minimize the bending energy, on surfaces of constant curvature, is considered in 
Refs. Il^, U, while interconnections with knot theory are explored in [|l2|. Adding the constraint on 
the enclosed area, as we propose here, provides a particularly fruitful generalization of this work. 
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Hamiltonians depending on extrinsic curvature have also been studied in statistical physics 1^, 
For space curves, an effective hamiltonian which depends on twist as well as bending provides a phe- 
nomenological description of stiff polymers, and in particular DNA (see e.g. Rcfs. Moving up one 

dimension, the leading term in the Helfrich hamiltonian, which describes the equilibrium configurations 
of lipid membranes, is proportional to the square of the traced extrinsic curvature K integrated over 
the surface, namely, the conformally invariant Willmore functional |ll|, |2^, |2^. In the variant of the 
Helfrich model proposed by Svetina and Zeks, the so-called bilayer couple model, constraints are placed 
on the surface, such as constant area, constant enclosed volume, and constant integrated mean curva- 
ture (the latter breaking the K —K symmetry of the problem) [^.0 However, the equations which 
determine the equilibrium are highly non-trivial higher-order PDEs. One motivation for introducing our 
variant of the elastica^ as a planar analogue of the above, is its potential as a toy model for understanding 
these higher dimensional membranes. At the simplest level, if a closed membrane possesses a symmetry 
along the z-axis, its profile at any z will be a fixed loop. The dimensionally reduced Helfrich hamiltonian 
is then the bending energy of this loop. It is possible to examine the loop hamiltonian exactly. While 
a lot is known about axially symmetrical closed configurations there is very little non-perturbative 
knowledge of the equilibria which exist when axial symmetry is broken: the loop configurations provide 
an analytic point of entry. 

An additional application of the planar problem was pointed out recently by Willmore: any closed 
loop can be exploited to generate an axially symmetric toroidal geometry. The equilibria of the loop 
hamiltonian on a surface of constant negative curvature can be mapped into equilibria of the Helfrich 
functional |p7|| . 

The feature of our model we will exploit is its integrability, very well disguised when the problem is 
cast with respect to the variables embedding the loop in the plane. Indeed, the determination of the 
curvature at equilibrium can be reduced to the study of the motion of a particle in a one-dimensional 
quartic potential. The curvature can be solved as a quadrature in terms of elliptic integrals. Remarkably, 
we discover that the problem possesses a second, far less-obvious, level of integrability: the embedding 
can be expressed as a local function of the curvature without any need of the two integrations one would 
have expected. As a consequence, the equilibrium configuration is given by a geometrical construction of 
Pythagorean simplicity. 

The condition of closure of the loop results in a discrete spectrum for a given constrained perimeter 
and area. We show that scaling can be exploited to fix one of these, the perimeter say. One branch of 
the spectrum then consists of configurations with an rt-fold symmetry (n = 2, 3, . . .) which inflate into 
a circle as the enclosed area is increased. If the area is evacuated these configurations will eventually 
self-intersect, tending to a limiting configuration of an inverted (i.e., negative area) circle, winding n — 1 
times and decorated with n infinitesimally small circles. The latter dominate, in this limit, the energy 
and give rise to a pole in the energy versus area diagram, the residue of which we evaluate exactly for 
each n. 

The remaining branches of the spectrum correspond to rather complicated self- intersecting pairs which 
bifurcate from a limiting configuration. While self-intersecting configurations are undoubtedly of interest 
mathematically, they are not of primary relevance to the physical problem we pose — we determine 
sufficient conditions for avoiding them (nevertheless, we describe a physical situation where they may 
arise) . 

In general, one would expect equilibria connecting perturbations to be governed by a fourth order 
differential equation. We show, however, that this equation can be reduced (by three orders) to a 
first order one. The latter is singular at the circular configuration, where bifurcation occurs. Analytical 
approximations for the energy as a function of the area are derived, in the limit of sinusoidal perturbations 
of the circle of frequency n. Combined with the pole mentioned above, these expressions provide a reliable 
sketch of the energy - area relationship for each n. 

^ We note that there is no genuine one-dimensional analog of its refinement, the so-called area difli'erence model . 
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The paper is organized as follows: Sect. || introduces the model and the shape equation that de- 
termines equilibria. We show how it can be reduced to the motion of a fictitious classical particle in 
a quartic potential. We address the behaviour under scaling of the shape equation and we show that 
the problem has essentially one free parameter. Moreover we analyze the loop statics; this provides a 
physical interpretation for the shape equation itself. In Sect. ^ we study the configuration space for this 
model. In particular, we describe the angle Oq by which the normal to the loop gets rotated in a full 
oscillation of the fictitious particle in the potential. Equilibria connecting perturbations is the subject 
of Sect. ^. These are determined by a complicated fourth order differential equation, which is used to 
obtain the purely geometrical construction mentioned above. In turn, the latter permits the reduction of 
the order of the original equation from four to one, as well as the derivation of a sufficient condition for 
non-self-intersections. Finally, we obtain analytical expressions for the energy as a function of the area, 
in some limiting regimes, which allow for a reliable sketch of its behaviour. In the appendices we collect 
various expressions, useful in the calculation of variations, we derive a recursion relation for the average 
of the powers of K and comment on an interesting Legendre transform in the space of parameters of the 
model. 

2 First Considerations 
2.1 Energy functional 

A closed loop, parametrized by s G [0, 1], is described by the embedding in the plane, 

x = X{s). 

The arc length £ along the loop is given by 

1/2 

ds' ds 







The most general expression for the energy of the configuration X{s) which i) does not depend on the 
parametrization, ii) involves no higher than two derivatives of X and Hi) is quadratic in these derivatives, 
is given by 

F[X] =a I d£K^ , (2) 



where K is the geodesic curvature, equal to the inverse of the radius of curvature at each point of the loop. 
We will take henceforth a to be equal to 1 — a is dimensionful, unlike its two-dimensional analogue, 
and this means we are measuring length in units of a. Let t be the unit tangent to the loop (transversed 
counterclockwise) and n its outwards normal. Then we have the Frenet-Serret equations for a plane 
curve^ (the prime denotes a derivative with respect to arc length £), 

n'^Kt, i'^-Kh. (3) 

In terms of the angle Q that fi makes with, say, the x-axis, K = Q' . 

To implement the constraints of fixed length and enclosed area, we introduce the constrained func- 
tional: 



F^[X] = F[X] +f^(^J d£-L^+a (^J^^ d^x ~ . 



(4) 



^Notice that the opposite sign convention for K is quite common in the literature. 
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fi appears as a lagrange multiplier enforcing the constraint fixing the length of the loop to some value, L. 
In the same way, —a is associated with the constraint fixing the enclosed area to the value A. We note 
that, in general, 

dF,{L,A) dF,{L,A) 

df—' " = 

In particular, if the area constraint is relaxed so that cr = 0, then dFc{L, A)/dA — 0. Looking at (||), one 
might be led to identify /i, a with the tension and differential pressure on the loop but this is only half 
true — we discuss the physical meaning of these parameters in Sec. We have not included, in our 
energy functional, a term proportional to the integral of K. This is because such a term is simply the 
total angle by which h gets rotated in transversing the loop once, equal to 2t: for a non-self-intersecting 
loop. 

We emphasize that the bending energy is not scale invariant (notice that neither is the dimensionally 
reduced Helfrich hamiltonian). An unconstrained circular loop will expand without limit. The natural 
scale invariant expression, Jd£\K\, though not in itself a topological invariant, does have a vanishing 
Euler-Lagrange derivative almost everywhere. Taken as an energy functional, its constrained configura- 
tions will be arcs of a circle (solutions of fiK = a ) joined together at curvature discontinuities in such 
a way as to mimimize the corresponding constrained energy function. This model will be considered 
elsewhere. 



2.2 Shape equation 

The Euler-Lagrange equations follow from extremizing Fc [X] , 



The reparametrization invariance of i^c[A] implies that the tangential projection of these equations is 
vacuous. The normal projection gives 




(7) 



(see the appendix for some relevant formulae). We note that with the identification, /i fi + Kq, Eq. (0) 
coincides with Eq. (7) of (see also [0). Is is worth pointing out that a loop on a surface of constant 
Gaussian curvature G satisfies the above equation with fi ^ fi — 2G. 

A plane curve is determined, up to rigid motions, by its curvature Q]. As is well known, our one- 
dimensional system is completely integrable. This is obvious because we have cast derivatives with respect 
to arclength. With respect to an arbitrary parametrization of the loop, first derivative terms associated 
with the one-dimensional Laplacian appear which conceal this fact. Writing (^ in the form 



^ ' 8 4 2 ' 



(8) 



we map our problem onto the problem of determining the motion of a fictitious particle in a quartic 
potential, K being the displacement of the particle and £ playing the role of time. The total energy E of 
the particle is conserved, 

^ = 0, F^lK-^ViK), (9) 

a fact that permits the expression of K' in terms of K and hence, of the arclength along the loop as an 
integral over K (via d£ = dK/K') 

. (10) 

^2iE-V{K)) 
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We emphasize the difference between the energy E of the fictitious particle, on the one hand, and the 
bending energy F of the loop, on the other. In particular, the configuration of least 'energy' in the 
analogue is not the minimum of the hamiltonian. 

The motion of the particle in the potential V is periodic. Closed loop configurations consist of an 
integer number of identical segments, each one corresponding to one full oscillation of the particle, and 
hence made up itself of two symmetric halves. Then, the condition of closure of the loop can be expressed 
as 

60 = —, 71 = 2,3,... (11) 
n 

where Qq is the angle by which the normal n gets rotated in one full oscillation of the particle, given by 

60 = 2 / = (12) 

(Kmin, Kmax denote the turning points). The equilibrium configurations have n-fold symmetry and a 
well-defined center. The value n = 1 is omitted in (|ll]) because it is special - see below. We examine in 
detail the resulting configurations in Sect. ||. 

2.3 Scaling 

Let us examine the behavior of Eq. (|^) under scaling of the position vector, X XX. We find 

Let X{£) correspond to some given solution of Eq. (^, with parameter values (cri,/ii). Then, a scaled 
solution with the same shape, XX{X£) and with (Li, Ai, Fi, Ei) (Aii, X'^Ai, X^^Fi, X^^Ei), is 
obtained by rescaling the multipliers as follows: 

1 1 
^ = Ma , (Ti ^ ^CTi = fJA ■ (14) 

Eliminating A we find the orbits of scaling in the (cr, /i)-plane 

^^ = — '^A • (15) 

Furthermore, an inversion in the origin, — > —X^^ can be identified with a rescaling by A = — 1, which 
maps a solution with a given a into one with —a. It follows that, for the purpose of identifying distinct 
configurations at least, one can set, e.g. a — 1 and scan the essentially one-dimensional configuration 
space by varying /i. For each value of /i, which fixes the form of the potential V{K), one has still to vary 
the energy E of the particle in the well so as to satisfy (|ll|), resulting, in general, in an infinite discrete 
spectrum "above" the point (1,/i) in the (cr, /i)-plane. 



2.4 Windings 

The closed configurations we wish to classify possess two topological invariants, both given as the winding 
numbers of maps — + S^. The first such map is the normal map of the loop, associating with each value 
of the parameter s (which ranges over the first above) the corresponding value of n (which ranges 
over the unit circle, the second above). For each rotation of s, n will generally rotate m times, with 
m an integer. Configurations with distinct m's are topologically isolated, the physical implementation 
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Figure 1: On the left: the potential V{K), for the values of /i, cr shown. Barred quantities, in this and 

— 0/2 

subsequent figures, are measured in units of ' , q being their length dimension. Also shown, with horizontal 
lines, are the energies that give rise to the n ~ 2,3,4,5 closed configurations, equal to 0.253, 0.278, 
0.298, 0.320 respectively (in units of /Iq and measured from the bottom of the well — compare with the 
p, = 0.5 curve in Fig. 0). On the right: the closed trajectories in the phase plane K — K' the ficticious 
particle, corresponding to oscillations in the potential V with energies E^'^\ E^^^ — the n = 3,4 trajectories 
lie in between the two shown. 

being the infinite energy barrier associated with the move ,_Q_^ . Non- 

self-intersecting loops correspond to normal maps with index m = ±1 (the converse is not true). 

The second map is from the s-circle to the closed orbit (~ S^) traced by the phase point in the plane 
{K, K'). For a configuration corresponding to n complete oscillations of the particle, the phase point goes 
around the closed curve ^K'^ + V{K) = E n times. Fig. |l|. There is a finite energy barrier separating 
configurations with distinct n's which prevents transitions between them. 



2.5 Forces and torques 

What can we learn about the loop configurations by looking at the equilibrium of forces and torques on 
an infinitesimal loop segment? First we ask, why would a rigid loop want to have area less than the 
maximal allowed by its length? One setup that supplies an answer is to imagine that our loop is actually 
the cross-section of an infinite cylinder, the interior and exterior of which are filled with an incompressible 
fluid. We start by filling the cylinder to its maximum capacity, this gives a circular loop. We then take out 
some fluid from the interior while making sure that the external pressure is large enough so as to prevent 
the formation of bubbles in the interior. Notice that, because of its rigidity, the cylinder takes up some 
part of the exterior pressure and only transmits to its interior a fraction of it, it's the differential pressure 
that crumbles the cylinder walls. The second of (||) then points to the identification of — ct with this 
differential pressure, taking into account that the latter points inwards. One has to be careful in applying 
the same argument to ^. The first of (^ seems to suggest that ^ is the tension of the loop but this 
would presuppose the possibility of tangential deformations — there is nothing in our energy functional 
that tells us how much energy these cost and, indeed, we have already assumed that the loop cannot be 
stretched or compressed tangentially. The derivative in the first of (^ is computed by comparing distinct 
loops with infinitesimally differing lengths, it does not refer to the deformation of a single loop. With 
this in mind, we now turn to loop statics. 

We denote by T(^o) the total force from the segment of the loop with t < £q to the one with I > 
and by r(€o) the corresponding torque. The latter is equal to —2K — this follows from our normalization 
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-ft(e + d£) . 


i ^ -fn{i + di) 




K{e)dt 
















(ll) 




Figure 2: Balancing of (a) torques, (b) normal forces and (c) tangential forces on a segment of the loop. K 
and K' are sketched positive, the resulting directions of torques and forces are shown. For Tt this data is not 
enough, we have assumed additionally > /i — this puts the center of the loop somewhere towards the left 
of each picture. 



of the bending energy F. Then, balancing the torques on a segment extending from i to i + d£ we get 

T{e)^T{e + de) + Tn{£ + de)de^o => T'-r„ = o, (i6) 

where T„ is the normal force, the reference point was taken at £ and the sign conventions arc shown in 
Fig. ^ The equilibrium of normal forces on the segment gives 

Tn{£)-Tni£ + d£) + Tt{£ + d£)K{e)d£-ad£^0 -T^^ + KTt - a = , (17) 

while the tangential components give 

Tt{£)-Tt{£ + d£)-Tn{£ + d£)K{£)d£ = ^ T/ + Xr„ = . (18) 



Notice how, in (|17l ) and ( [L8D , the curvature K is responsible for the tangential force at £ + d£ contributing 
a normal component at £ and vice-versa. Solving the above system of equations we find 



T = -2K . 



T„ = -2K' , 



(19) 



where, in the third relation, the integration constant was fixed to the value fi by the requirement that one 
recover the differential equation for K, Eq. (Q). We see from (|9|) that — ^ is the tension of the loop at its 
inflection points, if any. Notice also how (Q) is obtained upon substitution of the last two of ( p^ ) in (17), 
thereby identifying the physical origin of each of the terms in the former. In a forthcoming publication, 
where we extend our considerations to the case of a loop in space, it is shown how the above expressions 
for the forces and torque follow, in a model independent way, from an application of Noether's theorem. 

Our submerged cylinder model for the loop leaves no room for self-intersecting configurations, could 
there be any use for these? Imagine the loop made of superconducting material, in the presence of a 
uniform magnetic field perpendicular to its plane. In the limit where the magnetic flux due to self- 
inductance is negligible compared to the one due to the external field, the area of the loop has to be 
constant to keep the flux constant. One can then adjust the area by changing the magnetic field and, for 
small enough areas, self-intersecting configurations will appear. 
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3 Configurations 

3.1 Qualitative remarks 

For the purposes of this section, we may set, as explained above, cr = 1 in the expression for the potential 
V{K). Its critical points are given by the zeros of its derivative, i.e., by the roots Ki, i — 1,2, 3, of 

- ^IK -1 = 0. (20) 



We find 



where 



Ki = -K+ + iVSK. 



m 
12 



m 



K2 

m 
12 



m 



108+ 108. 



Mo 



(21) 



(22) 



where fj,o = 3/2^/'^. Of these, is always real while the first two are real only for n > /ig. Notice that, 
for any /x, the effect of the quadratic and the linear term in the potential is relatively important only in 
a neighborhood of the origin. For motions of the particle with sufficiently high energy E, the time spent 
by the particle in this region is negligible and the quartic term dominates the motion. As a result, we 
may conclude that Qq approaches zero with increasing E, due to the approximate symmetry K — > —K 
of the motion. The term linear in K spoils this symmetry and, for a > (as taken in this section), makes 
6o slightly positive for large E. As 6o descends, with increasing E, from this positive value to 0, it will 
cross all critical values 2n/n for n greater than some no. We expect therefore to encounter configurations 
with arbitrarily high n, for all /x — see Fig. Q 

On the other extreme, when the energy E is only slightly 
greater than a local minimum of the potential (at Ki or 
K3), and the particle oscillates around that minimum, we 
may approximate V{K) by a quadratic expression in K and 
find for 60 the limiting value 

V2 



e 



{max) 



2tt 



1,3, 



(23) 



where the superscript (max) is used because, as we shall see 
shortly, this is actually a maximum of Ooi(/i, E), for fixed /i. 
Starting from (pT|), one infers the limiting values 



-3__---2 M 


k 1 
— 1 

1 i = 3 
-0.5 , 

2 




1 *' 

1 m/Mo 

-0.5 

-1.5 / 
-2 ' 



©{max I / 
01 (m 



as well as 



e, 



{max) 
03 



(m 



0, 



(m 



= -27r, 

(24) 



Figure 3: Q'Si'^'^ip)' « = 1,3 



e 



03 



^(/i = 0) 



HI 



0{max) ( 
03 (M 



2ti. 



(25) 



A plot of Oq™°'^''(/x) is given in Fig. ^. We notice that, for /i > /io, there are no non-self-intersecting 
configurations corresponding to oscillations in the left well. 

We explain now why the value n = 1 was omitted in (pj]). The reasoning behind that relation was that 
a closed configuration should correspond to n complete oscillations of the particle. This is not necessarily 
so in the case of the circle. Any number of "oscillations" (of zero amplitude) will fit into a circle, including 
irrational numbers, and indeed, (^3|) shows that the circle corresponding to the particle resting at Ki, 

i = 1,3, is made of ^3 — -^I\f2 complete oscillations. 

To examine what happens for values of E that render the quadratic and linear terms important, we 
analyze the cases /i < ^0 and fJ. > Ho separately. 
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Figure 4: Qo3{E), E = (E — V (K'i)) j [i^, for various values of /2 = ^/ fxa- Also shown, with horizontal dashed 
lines, are some of the values of 9o3 that give rise to closed configurations. The corresponding values of can 
be read off as in the case /2 = 0.5, n = 2 shown — the (self-intersecting) configuration itself appears in Fig. ||. 
Notice how the Jl = —1, —2.5 curves miss the n = 2, rt = 2, 3, 4 configurations respectively. 



3.2 n < Ho 

The potential possesses only one minimum, at > given by ( |2l| ) (see Fig. |l]). O03 starts at Oo™"^"*, for 
E = V{K3), and for fi far from ^q, decreases monotonically to zero with increasing E. The configurations 
that do not appear are those for which 2 < n < 27r/8g'"''^'' as well as all with negative n. This condition 
determines a set {iJ,n\n = 2, 3, . . .} of critical values of ^, such that, for fi < fin, all configurations with 
n' ~ 2,3, . . . ,n are absent (the circle is, of course, always there). As fi approaches fiQ, the left wall of the 
potential develops a plateau that tends to the horizontal, as ^ ^ fiQ. The particle spends a relatively 
long time in this (negative K) region, which results in a negative bump in the Qq curve. As — > /io, the 
minimum of this bump tends to —00 — negative n configurations appear accordingly. Direct evaluation 
of the integral ( p^ ) gives ||] 



Qo3{fi,E) 



8g{aB + bA) 
A-B 



a2/c(fc) + |n(— 



,fc) 



(26) 



where fc^ = 



A^ = {a-hxf^a\ 
a + 4 



B 



A-B 
A + B 

(b-hf + al 

Kmax[iJ',E) 



0-2 

bi 
b 



_ bA-aB 

~ aB + bA 

= -i(a + 5) (27) 

= Kmini^J;E) 
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and /C, n are the complete elliptic integrals of the first and third kind respectively, given by 



lC{k) 



Vl - sin^ 9 



^0 fl 



(1 - a2 sin^ 6)^1 - P sin^ei 



(28) 



K 



A plot of 903 



The result (J^) holds for /i > /ig as well, with the appropriate choice of K„ 
for various values of fi, appears in Fig. ^. 

A couple of remarks are in order at this point. Consider the point A in Fig. ^ which is the intersection 
of the /2 = 0.5 curve with the n ^ 2 line, corresponding to an 7i ~ 2 configuration. Imagine now that p, is 
diminished continuously — the corresponding curve will move more or less downwards forcing A to move 
to the left. The corresponding E then diminishes which means that the fictitious particle oscillates in the 
potential well with smaller amplitude. When A hits the Ooa-axis, E is zero, the particle sits at the bottom 
of the well and the corresponding configuration becomes a circle. In other words, all configurations like 
the one corresponding to the point A can be continuously deformed to a circle. 

This is not the case though with the points Bi, B2 
in the same figure. Consider Bi — it corresponds to an 
n = —4 configuration for p, = 1.01. Imagine now that 
p, is diminished in a continuous way. As soon as it be- 
comes smaller then 1, the infinite negative pole of the 
corresponding curve is softened to a negative minimum 
which, for p sufficiently close to 1, still intersects the 
n = —4 line. Given that the /z-curves in Fig. |^ become 
almost horizontal at high E, we give an exagerrated 
sketch of the situation in Fig. ||. As the value of p is 
lowered, this negative minimum rises and, for a critical 
value of p, will just touch the n = —4 line, i.e., Bi 
and B2 collapse to a single point B. The configuration 
that corresponds to this point is not a circle, since E is 
positive and the particle oscillates with a finite ampli- 
tude, i.e., K is not constant. We conclude that con- 
figurations corresponding to points like Bi, B2 cannot 
be continuously deformed to a circle. Moreover, there 
exist bifurcation points, like B above, distinct from the 
circle. When p is lowered even more, the minimum of 

the /I-curve becomes positive and remarks similar to the above can be made about its points of intersection 
B'l, B'2 with sufficiently high n lines (see Fig. ||). 




Figure 5: Bifurcations: as p is reduced, the min- 
imum of the curve shifts upwards and Bi, B2 col- 
lapse to a point. Further decrease of p forces B[, 
B2 to collapse to a point. The corresponding con- 
figurations cannot be continuously deformed to a 
circle. 



3.3 M ^ Mo 

The potential posseses two local minima, at i^i < and K3 > (with V{Ki) > ^(iCg)) and one local 
maximum at K2 < (see Fig. |^). Depending on its energy and where it is started from, the particle is 
confined in the left well, the right well or visits both during every oscillation. 

We note the following, regarding the asymptotic behaviour of Ooi as ^ 00. In this regime, the 
linear term becomes negligible in V{K) and the integral giving 6oi can be easily seen to reduce (with 

—> x) to an integral giving (half) the period of a harmonic oscillator in a potential x'^ /A — ^x/2, which 
is independent of its energy E, as well as the 'constant force' /x. Another way to see the /j,-independence 
is by observing that, with tr = (which is equivalent to sending 11 to infinity), one lies on the ^-axis in 
the (cr, /z)-plane, which is an orbit under scaling, hence changes in fj, in this regime leave @oi invariant. 
Another novelty in this case is that, as the energy approaches V{K2) (either from above or below), the 
period of the motion tends to infinity, most of which the particle spends at K2. Assume we start the 
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particle at the turning point on the right with zero velocity and with just the right energy to reach the 
central maximum of V, E — V{K2). The corresponding curve starts in a counterclockwise sense with 
decreasing curvature, passes through an inflection point and acquires inflnite length spiralling forever 
clockwise as it approaches asymptotically a circle of radius I/K2, see Fig. ^ 

803 tends accordingly to —00, this corresponds to the negative poles in Fig. ^. As mentioned already, 
for V{Ki) < E < V{K2), So will have two branches, 801 and 803, corresponding to the motions 
confined to either well. Both of these branches will tend to —00 as i? ^ V{K2)- For E > V{K2) 
there is only one branch, which starts from —00 and approaches asymptotically zero (as E 00), after 
having reached a positive local maximum. A three-dimensional plot of 803 (/i,£^) is given in Fig.0 We 
present a representative collection of configurations in Fig. |8[ We also show, for reference purposes, some 
self-intersecting configurations in Fig. ^. 



4 Perturbations 



4.1 Equilibrium-connecting deformations 

We study here the following problem: given an equilibrium configuration X, find an infinitesimal defor- 
mation 5eX = en such that X + S^X describes a nearby equilibrium. Notice that S^X is not necessarily 
the displacement vector of some point of the loop — the latter will generally have a tangential compo- 
nent as well. Our motivation is to obtain new solutions from known ones, with the closure condition 
automatically satisfied. 

Using some of the formulae listed in the appendix, the linearization of the equilibrium condition (Q) 
gives that the deformation e must satisfy 

DEi{e) EE 2e"" + {5K'^ - fj:)e" + lQKK'e' + {12E -^K'^ + 6fiK^ + 10aK)e 

= -6,fiK - 5,a , (29) 

where the deformed configuration satisfies (0), with ji ^ ^ + S^^i, a ^ a + S^a. The energy will change 
as well, E ^ E + S^E. This appears explicitly in the linearization of (||), which gives 

i:>£;3(e) = 2K'e"' + {K^ - fiK - a)e" + 2K^K'e + K{12E-^K* + 2^K^ + 5aK)e 

= - ^K^6,n - K6,a - S,E . (30) 
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Figure 7: The angle Qo3{p,,E), for 0.5 < /I < 1.3, in steps of 0.05 and < i? < -8, in variable step. The 
/i = 1 curve has been shifted to p. = 1.01 in order to avoid the ambiguity in the number of roots of V. The 
negative pole in the curves for /I > 1 occurs at E = V{K2)- 

4.2 The cr-identities 

We introduce a number of functions on the loop, relevant in the study of the solutions of (p9|). Any vector 
field a defines the following two functions on the loop 

Qn = a ■ h , at = a ■ t . (31) 

The position vector X, in particular, supplies 

h = X-n, p=X-i, (32) 

both functions evidently dependent on the choice of origin. One finds 

a'„ = ~Kat , a[ = Kun , h' ^ Kp, p' = I - Kh , (33) 

where a is henceforth assumed constant. We will say that a function / on the loop is the generator 
of a transformation iff, under the latter, 6X ■ n — f. With this definition, a„, at are the generators 
of translations along a and perpendicularly to it respectively, p and h generate rotations and dilations 
respectively (both w.r.t. the origin). One easily verifies that DE/^{an) = DEi{at) = DE^{p) = 



while DEi{h) = 2^K + 3cr, in accordance with the scaling behavior of /i, a found earlier, Eq. (14). 
Experimenting a little with (^9|) we find, not without some surprise, two more solutions 

DEi{K') = , DE^ {<7'\K^ - m)) = 2^*/^ + 3(7 . (34) 

Noting, additionally, that when the origin is at the center of the loop, the zeros of p coincide with the 
extrema of K, as well as the coincidence of the extrema of h and K^, we make the ansatz 

h^a-\K^-n) + f, p^2a-'K' + g, (35) 
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Figure 8: Non-self-intersecting configurations for /i = —0.5, a — 1 and n — 2,3,4,5. Also shown, for each 
configuration, is the circle of radius Xq (see Eq. (psj)), passing through its inflection points. 



with /, g to be determined. Substituting in (p3|), we find /' = Kg, g' = —Kf, with solution 

/ = a cos(e - ^o) , g = a sin(e ~ (j)o) , (36) 

where a, (f>Q are arbitrary constants. One recognizes /, g to be, respectively, the functions a„, at corre- 



sponding to the position vector a of the center of the loop, then (35) state 




(37) 



a remarkable geometrical property, the implications of which will occupy us shortly. Notice that neither 
of (^7|) survives in the cr — > limit. Unless otherwise stated, we will take a equal to zero in what follows. 

One might wonder whether (|28| ) admits other polynomial solutions in K, apart from the second 
of (^J). To investigate this, we rewrite ( p9|) in terms of derivatives w.r.t. K, using ^ = K' and find 
for DE{K'^) the leading term {^n^ + jn? + — 3n — |)if"+'*, with only positive integer root n = 2. 
We conclude that no other polynomial solutions of (k9) exist. 
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4.3 Geometries 

We are now in a position to give a purely geometrical construction of the equilibria. Indeed, starting 
from = + and using (|3^) , we find 



-Xi. = 4(7-^ K, 



Xo = a-^y/8E + f^, 



(38) 



where, in the derivation, use was made of the first integral, Eq. (j^). This remarkable formula expresses 
the embedding completely in terms of K{£). We emphasize that the shape is obtained directly from K, 
given {a,fj,,E), without any integration. Note that both the relations ( |3^ ) follow by taking derivatives 
with respect to arclength of ( p8| ) . It follows from (^8|) that the moment of inertia J of the loop around 
an axis perpendicular to its plane and passing through its center, is given by 



I = LX^ + 8n<T-^ 



(39) 
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Volume = 4cr 




p = 



l_ 

K 




Figure 10: Geometrical construction of the equilibrium curve. The parallelepiped has square base with side 
AP and volume 4a^^. Its height p gives the radius of curvature at P. The part of the curve lying in the 
interior of the circle is constructed similarly. 



To construct geometrically the equilibrium curve (for given a, /i, E), we draw a circle with radius Xq 
and from a point P outside it, bring the tangent PA to the circle (see Fig. |l^). Taking as base the square 
of this tangent, we construct a parallellepiped with volume 4cr~^, the height of which is the radius of 
curvature of the equilibrium curve at P — this follows from (|3^), written in the form {X^ — X§)p = 4cr~^, 
where p is the radius of curvature at P. To achieve closure of the resulting curve, one has to start with 
a particular slope, given by either of (p^). In this way, one constructs the part of the curve lying outside 
the circle (the latter intersects the curve at its inflection points, if any) — the interior part, present only 
if the ficticious particle reaches into the negative K region, is constructed in a similar manner. 

Referring back to our expression for the force T{£), Eq. (19), we realize that ( |37|) implies that f{£) is 
orthogonal to X{£) while T{e) = aX{£), i.e., 



T{£) ^a{hxi)x X{e), 



Tf — ah . 



(40) 



with T{£) the force on the part of the loop pointed to by t. Notice that f is generally compressing but 
changes nevertheless to true tension at the points were X is tangential to the loop. Also, p has acquired 
a direct physical interpretation as a consequence of (^): apd£ is just the torque, w.r.t. the origin, of the 
force due to the pressure on a segment of the loop with length d£. Then the vanishing of the total torque 
on the loop is guaranteed by the fact that p is a derivative. In fact, one may derive a compact and rather 
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pleasing formula for the torque due to pressure on any segment of the loop, like the one defined by A, B 
in the sketch of an rt = 4 configuration of Fig. 11. 

Suppose we take A as reference point, then (37) is valid with a connecting 
A with the center O of the loop. We get 



'AB 



.e., 



a dip ^ a I d£{2a-^K' + at) = 2{KB-KA) + aa- 
Ja J a J a 



de 




(41) 



Aj 

















Moving the reference point simply moves one endpoint of a. For a — 

JO) 
'AB 



Figure 11: Computing the 
torque due to pressure on 
the arc AB 



(torque w.r.t. O), t^b = 2(Xb 
4.4 Self-intersections 

Given that some potential applications of our model exclude self- intersecting configurations, we look 
now for sufhcicnt conditions for non-self-intersection. We restrict our attention to self-intersections that 
can be reached by a continuous deformation of non-self-intersecting configurations. In other words, we 
consider a one-parameter family of configurations X{(.^t), t G [0, 1], continuous in i, such that, for every 
t, the corresponding curvature satisfies (|^) and we take X{l,Q), X{t,\) to be non-self-intersecting and 
self-intersecting respectively. Then, we observe that as t varies from to 1, one necessarily encounters a 

"kiss" : The position vector Xa is along the axis of symmetry of the lobe and tangential to the loop 

at A. The third of ( ^0|) then shows that the force at A is normal to the loop, while the second of (|40| ) 
gives its magnitude as 

Ta - (jXa ■ (42) 

The tangency of X at A implies that /i^ = and hence, using ( ^ ) once more, we get (assuming that K 
is negative at a kissQ 



Ka = -Vm, 



(43) 



for all kisses, regardless of the order of the configuration. We conclude that a sufficient condition for 



non-self-intersection (of the type defined above), is fi < 0. (|43| ) gives K'^ in terms of cr, ii, E 



K'\^2{E+\^?~\a^ 



On the other hand, from (|3^) we find that 

X\ = cr-\8E + n^)-A(j-^^. 



(44) 



(45) 



The tangency of Xa though means that 2a ^K'^ — pA — Xa and the above two equations then give 
(setting a — 1) 

£^' + ^A''-4^/M-^)i?+^(M'-4VM)(/i'-4VM-l) = 0, (46) 



This assumption is true for all configurations we have studied numerically. 
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with roots 



1 ^ 1 n 1 

^^-2^-8^+8 



(47) 



For the particular case n = 2, K'j^ — and we get E — £2- We have seen in Sec. ^ that the relation 80 = 
27r/n defines a curve in the (fi, E) -plane, consisting of all parameter pairs giving rise to a configuration 
of order n. The intersection of that curve with the ones just written above, consists of the points (/i, E) 
giving rise to a kissing configuration of order n. 

Consider how (^o|), (|4l| ) guarantee equilibrium in some particular examples. First, look at an even-n 
configuration, say, n = 2: A^^^^^^ B The total force due to the pressure, pushing together the 

two halves in the sketch, is ct|AB|. ( ^o|) says that the tension at A, B (purely tangential, compressing) 
is a\AB\/2, thus leaving each half at rest. As a second example, consider the lobe defined by a self- 
intersection: The total force due to the pressure on the lobe is zero, and ( ^o|) says that the forces 

from the rest of the loop, on the two ends of the lobe that meet at A, are opposite (with direction so 
as to keep the lobe closed). As a further check on our results, one can verify the balancing of torques 
using (|4|) on, say, the right half of a kiss. 

4.5 Connecting equilibria 



(48) 

(49) 
(50) 
(51) 
(52) 



There is more to be derived from (|3^). Taking 5^ on both sides, we find 

DE[°'\e) = 2e" + {iK'^ - = i^^^Xl - 2K(j-^5^a , 

a considerable improvement over (|2^)Q. Moreover, by varying the first of (^^ we also find 
DE;^ ^(e) = 2Ke" - 2K'e' + {a + 2K^)e = -a-\K^ - fi)S,a - S^n . 



Comparison with (48) leads to a first order equation for e 

DEi{e) = AK'e' - AK"e = caA'^ + ciK + cq , 
where the constants Ci are given by 

C2 = — 2CT^^5e(T , Ci = (tS^Xq , Cq = 2(5e/X — 2/XO'^^(5£(T . 

Noting that ^ — (A')^^^, ( |50| ) can be written in the form 

8(A - V)e + We = C2 A^ + ci A + cq , 



where the dot denotes differentiation w.r.t. A and V — V{K). This latter equation is readily integrated 
to give 



(53) 



One may add an arbitrary amount of A' (rotation) to this (but not a„ or at, since these move the origin). 
To get particular solutions from (^3|), wc need to specify the direction of the deformation in the (cr, /i)- 
plane, i.e., the ratio b^ytji^a. It will prove convenient for our further analysis of (^3|), to rcparameterize 
the (cr, /z)-plane introducing new coordinates (A, f ) via 




A(g, y)^a ' , i[a,y)^a '-fi. 

^Notice, however, that is only valid when the center of the loop is at the origin. 



(54) 
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A point P with coordinates (cr, ^) 7^ (0, 0) lies on a unique scaling orbit, which can be specified by the 
/i-coordinate of its point of intersection P' with the a — 1 line, this is the value of ^ for P. One can get 
now from P' to P by scaling by A (see Fig. |l^). The obvious advantage of these new coordinates is that 
the scaling orbits are constant-^ lines. The dependence on A of an arbitrary quantity iS'(A, ^), with length 
dimension q, is ^(A,^) — A'^S'(^), where we denote by a tilde the remaining function of ^. It follows that 
^ = qX^^S so that S^S = q\~^SStX + X'^S'S^^, the prime denoting here differentiation w.r.t. ^. 

We verify that ( |53| ) gives e ~ /i for scaling. Eq. (pT]), written in 
terms of (A,^), gives 



C2 
Cl 

Co 



= A-2 



( WE 



2^' (5,A + A-M8^' + 2CUeC (55) 



= 2X-^^S,X + 2X-'S,^, 



where E{X,^) = X~'^E{£_) is the energy that guarantees closure of 
some particular configuration. An increment S^X corresponds to 
scaling by a factor 1 + S^X/X, hence the corresponding e should be 
e = {S^X/X)h. Putting S^^ = in ( [56| ) one determines the Ci for pure 
scaling, then computing d/dK{{K^ — ii)/K') one finds that the 
and terms in the numerator cancel and one recovers the integrand 
in ( |53| ) with just the right c's. 

Integrating the first of (p7|) w.r.t. M, and using 



A 










Xp' I 




— 1 — 1 — > 



(7=1 a 



Figure 12: Definition of the 
coordinates A, ^. 



A=i 
2 



Mh, 



we get 



from which 



F = plL + 2a a , 



S^F = LSe^L + IJ-S^L + 2AS^a + 2a5eA 
follows. On the other hand, direct variation of F, Eq. (^, gives 

5^F = —fiS^L — aS^A . 

Comparing with (p8|) we find 



LS^fi + 2AS^a + 3aS^A + 2fiS^L = 
For a length-preserving deformation, ScL — 0, and ( |60| ) reduces to 

LS^fi + 2AS^(T + SaS^A = . 



0. 



(56) 

(57) 
(58) 

(59) 
(60) 

(61) 



We note in passing that the various differential operators we have defined above are related in the following 
way 

[DEiie)]' = 2K'DE2^''\e), 
K[DEi{e)]' = 2K'DE2^''\e)+K'DEi{e), 

DEsie) = E'lDE^^^'Ke)]' - E'DE^'-^He) -\{K^ - t^)DE,{e) , 

K'DEiie) = 2[Z?^3(e)]'. 

It follows that any perturbation that satisfies Eq. 



J) will necessarily satisfy all the higher order ones. 
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4.6 Perturbing the circle 

Referring to our cylinder model of the loop (see Sec. |2.5| ), we consider here the following thought exper- 
iment: we start with the loop "filled" to capacity, i.e., in a circular configuration, and start deflating it 
by removing fluid from its interior while keeping its length fixed. We would like to follow the evolution 
of its shape, for distinct n's, and derive the limiting form of the bending energy _F as a function of the 
area A near the circular extreme. Going beyond the cylinder model, we would also like to allow for 
self-intersections — what happens if one just keeps subtracting area? 
Consider the following perturbation to a circle of radius R 

p{e) ^ Ra{l + asm{n0)) , a<l. (62) 

The line element and curvature in these coordinates are 

di-^/^T^de, k=p100-, (63) 

where the dot here denotes differentiation w.r.t. 6. Substituting (|6^ ) in the first of (|63|), integrating over 
9 and requiring the total length to be equal to the initial value 2ttR, we find 

R^ = R{l-!^)+0{a'), (64) 
while for the curvature the second of ( |63| ) gives 

K=^(l + a{n^ ~ 1) sin(?i6l)) + 0{c?) . (65) 
R 

We substitute the above expression in the differential equation for Eq. (|^), and demand that it be a 
solution to Oia). The constant and Oia) terms respectively give 

1 - i?V - i?V = , 271^ -I- (i?V - 5)n2 -I- 3 - = . (66) 

Notice that the first relation is simply the statement y(l/i?) = while the second also follows from 
with e as in (|6^). For this perturbation the differential equation cannot be satisfied to 0{o?). From (6^) 
we get 

^=;^("'-l), ii=^{i~2T?). (67) 

During the deflating process, both a and /i will vary as functions of A. The point in the a-p, plane 
corresponding to the configuration will trace out an orbit, starting at the above points (for each n), all 
of which lie on the line i?'^(T + R}[l ~ 1. For F and A we find 

^(1 + 1(^2 _ 1)2^2^ ^0(^3)^ A = 7ri?2(l-i(n2„l)^2)+0(^3)^ (gg) 



Notice that, near the circle, 



bF 2 



n^-l) 



in agreement with ( |59| ) (for = 0). ks, A keeps diminishing, our numerical analysis shows that the 
configurations start self-intersecting, giving rise to regions of negative area — this scenario is sketched in 
Fig. Notice how the winding numbers (-1-1 for the little circles, —1 for the big one) add up correctly 
to match that of the circle to which this configuration can be continuously deformed. This observation 
points to a feature of the above m = 1 limiting configurations already alluded to in the introduction: the 
big circle in the configuration of order n actually winds around itself n — 1 times. For example, to move 
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from one lobe to the next along the loop, in the n = 3 configuration, one must travel an angle of 47r/3, 
not 27r/3 (a glance at Fig. ||.b reveals how this comes about). The limiting shape in this direction then is 
a circular one, with radius slightly less than R/ [n — 1) and n little circles attached to it. In this extreme, 
F{A) is dominated by the little circles and assumes the limiting form 

^A)^'^ ^, , (69) 

where the area is counted with the appropriate multiplicity due to the winding {e.g., plus three times 
the area of the little circle, minus twice that of the big one, in the n = 3 case). A sketch of F{A) that 
interpolates between (^) and (|9|) is given in Fig. 0. 

5 Conclusions 

In this paper, we have studied the equilibrium configurations of elastic planar loops with constant area 
and constant length. The condition of closure of the loop gives a discrete spectrum of configurations, lying 
along several branches in parameter space. We focus on non self-intersecting loops that inflate to a circle 
when the enclosed area is increased, due to their relevance as toy model for two-dimensional membranes. 
For this branch, starting from analytical expressions in the relevant limit cases, we obtain a reliable 
sketch of the dependence of the energy on the area. The equation that determines equilibria connecting 
perturbations led us, rather unexpectedly, to the a-identities (^7|), from which a novel geometrical con- 
struction of the (constrained) elastic loop followed. In turn, these identities permit the expression of an 
equilibrium connecting normal deformation in terms of the curvature and the appropriate perturbations 
of the parameters. 

In future work, we plan to address the important issue of the stability of the equilibrium configurations. 



In light of the complicated structure of the fourth order differential operator appearing in (29), this 
appears to be a non-trivial task in the general case. It is indicative of the intrinsic complexity of the 
question that even the relatively simple case of the figure eight configuration requires the elaborate 
analysis of p^ . 

Another interesting issue is the analysis of thermal fluctuations in this model, which would provide an 
analytical counterpart to the study of Ref. of two-dimensional vesicles, using Monte Carlo techniques. 
In particular, there is a similarity between some of the configurations of Fig. 2 of the above work and 
the configuration of Fig. ^.a in this paper. It is also interesting to note that the shapes of Fig. ^ closely 
resemble a top view of the starfish vesicles of Ref. which are almost planar. Finally, configurations 
similar to ours appear in where our problem has been approached from a functional analytic point 
of view and important existence results have been derived, as well as in [p| . 
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After the completion of this paper, we set out to study its most natural extension: geometric models 
for loops in space. We became aware of a large body of literature that explores the interconnections among 
hierarchies of functionals of the geometry of a curve, their corresponding generators of rigid motion and 
integrable systems, such as the KdV equation. We have found Ref. an excellent point of entry to the 
subject. We also realized that our (T-identities can be obtained by adapting the cylindrical coordinates 
used by Langer and Singer in the analysis of buckled rings. We expect that these explorations, apart 
from their intrinsic interest, will contribute to a deeper understanding of two-dimensional membranes. 
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A Some Useful Formulae 

Under an infinitesimal deformation of the loop along the normal, X —t X + 5^X = A + e{l)n^ we find 
= Kedi, 6 J = e'n, S,K = -e" - K^e 

5,4 = -Ke4-„, S,n = ~e'i, S,K' = -e'" - Kh' ~ SKK'e 

at at 

S,h = e-pe', 6eP = e'h, 5,K" = -e"" - K^e" - bKK' e' 

-{6E - + ^pK^ + 5crK)€ 

(we have used the e.o.m and the first integral, Eqs. (j^) and (||), to express K" and {K')^ in terms of K). 
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For the derivatives of /i, p, K we find 
h' = Kp 

h" = -K'^h + K'p + k 

h'" = -mK'h+{-^K^ + ^K + ^)p + 2K' 



•III! 



p' = I- Kh 
p" = -K'h - K^p 

P'" = ilK'-!lK-^)h-3KK'p-K^ 



nil 

p 



i^K^K' - ^K')h + (i^X"* - l-nK"^ - 5aK - 6E)p - 5KK' 



^ 8 4 2 ' 

K" = ^{-K'^ + fiK + a) 

K'" = -Ik^k' + ^K' 
2 2 

2 4^ 4 M ^ 4 

B Averages 

Consider the quantities Wm defined by 

W„= / K''di^2 / ^"i77. (70) 

Lq here is the length along the loop corresponding to one full oscillation of the particle (we will use Fq 
later with a similar meaning). Wn is then (proportional to) the average of along the loop. Starting 
from 

d I 
if™.,. dK 
one may derive the recursion relation 

-i(n + 2)W^„+3 + J(n + l)V7„+i + ^(2n+l)W„ + £;nW^„_i = 0, n = 0,l,2,..., (72) 

which permits, in principle, the calculation of the average of any power series in K , for a given equilibrium 
configuration. In particular, 

W3=nWi+aWo, W4 ^ ^W2 + 2aWi + ^EWo . (73) 

Referring back to Fig. |^, we notice that the area S enclosed by the orbit of the phase point, for some 
given (T, /i, E, is given by 

S{a, n,E)^2 / K'dK . (74) 
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Differentiating w.r.t. a and using the fact that the integrand vanishes at the endpoints, we find 

e„^2f. (75) 

Length and bending energy per full oscillation also follow by simple differentiation 

dS dS 

Notice that these quantities are defined whether or not the closure condition is satisfied, i.e., they are 



functions of the three independent variables cr, /i, E. Combining (|74| ) and (73) above, we get 



S^^EL+^-^^F+]^cJQ. (77) 

One may regard Lq, 8o and Fq as new coordinates in the space of (not necessarily closed) configurations 
— the change of coordinates is non-singular, except for some special points, and is given by a Legendre 



transform. Then the closure condition, Eq. (11), restricts to the configurations that lie on the Lq-Fq 
plane, with 6o = 27r/n. 
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